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PACS 83 . 50 . Lh - Slip boundary effects (interfacial and free surface flows) 

PACS 61.25.Hq- Macromolecular and polymer solutions; polymer melts; swelling 

PACS 82 . 20 . Wt - Computational modeling; simulation 

Abstract. - Using molecular simulations, we study the properties of a polymer brush in contact 
with an explicit solvent under Couette and Poiseuille flow. The solvent is comprised of chemically 
identical chains. We present evidence that individual, unentangled chains in the dense brush 
exhibit cyclic, tumbling motion and non-Gaussian fluctuations of the molecular orientations similar 
to the behaviour of isolated tethered chains in shear flow. The collective molecular motion gives 
rise to an inversion of hydrodynamic flow direction in the vicinity of the brush-coated surface. 
Utihsing Couette and Poiseuille flow, we investigate to what extend the effect of a brush-coated 
surface can be described by a Navier slip condition. 
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• Introduction. — Grafting polymers onto surfaces is 
a versatile and stable method for controlling wettability, 
lubrication, adhesion and surface interactions [1]. Brush 
coatings are utilised for biocompatibilisation or as an- 
tifouling coatings in microfluidic devices. They are em- 
ployed to stabilise colloids or, in form of interfacial layers 
of copolymers, to prevent coalescence of droplets in poly- 
iner blends. ^ Dense polymer brushes give rise to fasci- 
nating dynamical properties like a dramatic reduction of 
friction between two brushes sliding past eacli otlier [2] , in- 
completely understood collective dynamics [3], additional 
dissipation mechanisms of droplets moving on brushes due 
to the the deformability of the soft elastic surface [4] , and 
large effective slip lengths [5]. 

Simulations of dense polymer brushes with an explicit 
solvent [6-9] are computationally very demanding and pre- 
vious studies of brushes in shear flow [10-14] conceived the 
brush as a porous medium [15] utilising Brinkman's equa- 
tion [16] to quantify the penetration of the solvent flow into 
the brush. They observed that density profiles normal to 
the surface are not strongly affected by shear flow but the 
average molecular conformations are tilted towards and 
stretched along the direction of the flow. 
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^Copolymers at interfaces in blends can laterally move, but the 
ends of grafted chains in a brush are immobile. 



Isolated grafted chains in shear flow already exhibit 
intriguing dynamics [17-20]: Thermal fluctuations away 
from the surface expose the non-grafted chain end to a 
faster flow and lead to an elongation of the chain in the 
flow direction. The stretched chain rotates back towards 
the grafting surface and contracts. Related effects have 
also been reported for the chain ends in swollen brushes 
under strong shear and they have been attributed to the 
non- linear elasticity of the macromolecules [14]. One hall- 
mark of this dynamic behaviour of isolated chains are 
strongly non-Gaussian distribution functions of angular 
orientations [18,20]. Explicit hydrodynamic interactions 
[19, 20] appear not to be crucial for this cyclic tumbling 
motion. 

Using Single-Chain-in-Mean-Field (SCMF) simulations 
and Molecular Dynamics (MD) simulations of a brush 
in contact with an explicit solvent of identical chain 
molecules we demonstrate that the tumbling motion ob- 
served for isolated chains persist as one increases the graft- 
ing density and that the individual motion of the tethered 
molecules results in the inversion of the collective flow di- 
rection at the surface. The effect of the brush cannot be 
described by a Navier slip condition [21-23]. 

Model and simulation technique. Brush and 
melt polymers are represented by a bead-spring model 
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with N = 32 segments per molecule. Harmonic bond- 
ing forces Fb with spring constant 3{N — \)kBT/E^ act 
between neighbours along the chain. denotes the 

mean square end-to-end distance, fc^ Boltzmann's con- 
stant and T temperature, respectively. Repulsive in- 
teractions between segments are catered for by penalis- 
ing fluctuations of the local density by an energy [24], 
Hnb = ^^po J d^r ((/)b + 4>m- if where ^b and (jj^, 
denote the normalised density of brush and melt, and 
kN = 50 characterises the inverse compressibility. The 
segment number density, po, determines the invariant de- 
gree of polymerisation, = (poR^/N)'^ = 128^. We study 
the flow of a polymer melt over a brush in a slit pore of 
width Dx- The lateral extensions are Ly = L^, and pe- 
riodic boundary conditions are applied in y and z direc- 
tions. At a; = and D^, there are repulsive, structureless 
walls ^ that repel segments via a potential ^ of the form: 
Vw(Aa;) = Aw exp(-Aa;V2e^) with AwA^ = ee.gfcsT, 
€yf = O.lbRe, and Ax being the distance of a segment 
from the surface. The first segment of the brush is grafted 
at Axo = O.lRe- Couette flow is imposed by moving the 
grafting points of the brush at the top surface at fixed 
velocity in y-direction; Poiseuille flow is generated by ap- 
plying an external body force, Fex, to all segments in y- 
direction. 

In SCMF simulations [25,26] one considers a large en- 
semble of explicit chain configurations that independently 
evolve in an external field. This field mimics the effect of 
the non-bonded interactions W{r) = (5Hnb/'5[/5o0b(r)] + 
Kv. Its gradient gives rise to a force Fnb = — Viy. 
Chain segments evolve via Smart Monte Carlo (SMC) 
moves [27, 28] with trial displacements 

ArtHai = (C(v) + F) At/C + ^^2kBTAt/C (1) 

F = Fb-l-Fnb+Fox denotes the forces acting on a segment, 
C is the segmental friction,^ and ^ is a Gaussian random 
number with zero mean and unit variance. We set CA — 1. 
^ (v(r)) denotes the hydrodynamic velocity field. Chains 
exhibit Rouse-like dynamics characteristic for unentangled 
melts [29] and the Wcissenberg number is defined as Wi= 
■yd^NRl/ (Sn'^kBT) where 7 denotes the shear rate. 

The fluctuating, external field, W, is calculated after 
every time step, and the densities are constructed by as- 
signing the instantaneous particle positions to a three- 
dimensional grid with linear spacing AL — i?e/6 [26]. The 
hydrodynamic velocity, (v) , however, represents the aver- 

^The melt perfectly slips over the bare surface (without brush). 
At the considered grafting density, however, the velocity at the wall 
vanishes (microscopic stick, see Figs. 2 and 4). Thus, our results do 
not depend on the microscopic boundary condition (slip or stick). 

repulsive VW impedes the SMC to propose moves that would 
penetrate the hard walls and be rejected. 

*The segmental friction, ^, is not determined by the intermolec- 
ular forces like in a MD simulation but it is a parameter of the 
algorithm. 

^For Brownian dynamics the equations of motion are invariant 
under C C' = -^C ^■•id t — » t' = At and we measure time in units of 
CNRl/kBT 



age flow field and must not fiuctuate. First, we directly 
calculate an instantaneous velocity, v, = Avi/At, of seg- 
ment i from its explicit displacement [30] during a SMC 
step (in order to retain spatial resolution) and assign it to 
the grid. Then, we add the instantaneous velocities of all 
segments, average over a time period, T, and normalise 
by the local density to obtain the average local velocity 
(v(r)). This procedure ensure that the average force (F) 
vanishes. This time averaging procedure limits the sim- 
ulation technique to stationary or slowly varying flows. 
For the density utilised in our simulations a time interval 
T > 800CRl/{NkBT) is sufficient to eliminate fluctua- 
tions of the velocity field and to yield accurate predictions 
in equilibrium. 

The computational scheme is similar to dynamic self- 
consistent field techniques [31-33] or self-consistent Brow- 
nian dynamics simulations [12,14,34]. With the SMC algo- 
rithm we systematically approach the Brownian Dynamics 
limit (At ^0). In the present work, we use the time step 
At = OmCRl/{NkBT) which closely mimics Rouse dy- 
namics. The time step is chosen such that the self-diffusion 
coefficient and the viscosity in the bulk are close to the 
Rouse limit, = /cbT/(CA) and T]BRe/{CNVM) = 
1/36, respectively. The acceptance rate of SMC moves 
is about 83%. We have explicitly verified that a ten times 
smaller time step yields very similar results. Moreover, our 
simulations capture density fluctuations [26] and, most im- 
portantly, we do not approximate the effect of the brush 
on the solvent flow by that of a porous medium [15, 16] 
but use an explicit solvent instead. 

To validate the dynamic SCMF simulations, we ad- 
ditionally perform MD simulations of a standard bead- 
spring model [9, 35]. Grafted and melt chains com- 
prise A^ = 32 segments interacting via a purely repulsive 
Lennard- Jones potential with length and energy scales, 
ctlj and elj, and cut-off To = v^ctlj. Neighbouring seg- 
ments along a chain arc bonded together by a Finitely 
Extensible Non-linear Elastic (FENE) potential. The tem- 
perature, fesT/cLj = 1.68, is maintained by a DPD ther- 
mostat [36]. The MD simtilations correspond to a dense 
brush, aRl = l.Sv^/V", of short chains, J\f = 6^. 

Results. — In Fig. la wc; present the density profile 
across the slit pore for Dx — lORe and grafting density 
(Ti?g — VM for Couette and Poiseuille flow as obtained by- 
dynamic SCMF simulations. At this moderate grafting 
density there is a broad interface between brush and melt 
("wet brush") and, in agreement with previous studies [14], 
the profile normal to the surface does not depend on flow 
for small shear rates. The depletion at the wall stems from 
the repulsive segment-wall potential. 

Measuring the forces, F^raft^ ^j^g^^ g^^^ ^j^g grafted 
segments in Couette flow we obtain the shear stress 
UxyRl/kBT = 29.3 for Couette flow. Alternatively, we 
can estimate the intramolecular stress of the melt at 
the centre of the film via the mean-field approximation 
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Figure 1: a) Density profiles across a slit of width = 10-Re 
and lateral extent Ly = — 4ile for Poiseuille flow {F^^N — 
QmksT/Re, symbols) and Couette flow {■y^CNRl/ksT = 5, 
Wi= 0.247, lines), b) Kramer's stress for Poiseuille and Cou- 
ette flow. The horizontal dashed line marks the stress obtained 
from the forces on the grafted ends for Couette flow. 



(Kramer's formula [29] ^) a^yRl/ {kBTVAf) 



l)(0b + 0m) R'i /[3(N-i)] '^here and by denote the dis- 
tance between bonded segments perpendicular to the sur- 
faces and along the shear direction, respectively. This es- 
timate is shown in Fig. lb for Poiseuille and Couette flow. 
The stress in Couette geometry agrees well with the result 
obtained from the force on the grafted segments. From 
the stress at the centre of the film and the velocity gradi- 
ent (cf. Figs. 2 and 4) we can estimate the dimensionless 
viscosity, r]RJ{(NV^) = 0.030 and 0.03f for Poiseuille 
and Couette flow, respectively. 

Fig. 2 presents the velocity across the film for Poiseuille 
flow. The profile at the centre is parabolic, 



(x) 



2r] 



(x - xp) (Dj; - xp - x) (2) 



as expected from macroscopic hydrodynamics, xp denotes 
the position at which this parabolic profile extrapolates to 
zero. The flow penetrates into the brush [f5] and details 
of the flow in the vicinity of the surface are shown in the 
panel (a) of Fig. 2. Unexpectedly, the hydrodynamic ve- 
locity profile in the vicinity of the grafting surface becomes 
negative, i.e., the flow direction is opposite to the driving 
force Fox. This inversion of the collective flow direction at 
the grafting surface is a consequence of the individual mo- 
tion of grafted chains, i.e., the tumbling motion of grafted 
chains under shear [f4, f7]. 

To rationalise the observed inversion of flow at the 
surface, it is important to distinguish between the time- 



^We utilise tlie same spatial assignment of the stress as for the 
density but other schemes [e.g., J. H. Irving and J. G. Kirkwood, J. 
Chem. Phys., 18 (1950) 817.] that are locally more accurate can be 
employed. 




Figure 2: Mass averaged, total velocity profile, (vy), and veloc- 
ity of brush and melt segments for Poiseuille flow. Panel (a) 
highlights the total velocity at the surface and shows details of 
4'hVh,y for a system that additionally contains 1% tracer parti- 
cles (i.e., 6554) which have the same extension as a monomeric 
unit and friction coefficient C = 1- Note that the velocity pro- 
files of the tracers coincides with (v). Panel (b) shows the 
velocity profiles across half of the film. The dashed line is a fit 
to macroscopic predictions (cf. Eqs. (2)). The arrow marks the 
positions, xp, at which the extrapolation of the macroscopic 
hydrodynamic velocity profile vanishes. 
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Figure 3: Sketch of the molecular motion of the grafted 
molecules of the brush, the free chains of the melt and the 
smaller tracer particles. On the left, brush molecules in two 
different stages of their tumbling motion are depicted and the 
velocity of the end-segment is highlighted by an arrow. This 
graph illustrates that the collective velocity of the brush seg- 
ments at a specific distance from the surface need not to vanish 
although the time-averaged velocity of each segment of grafted 
chains is zero. In the middle, the forces that act on the tracer 
particle and a segment of a free chain near the surface are illus- 
trated. On the right hand side, the collective, hydrodynamic 
velocity profile, (v), of all segments is shown. 



averaged velocity of individual segments and the position- 
dependent, hydrodynamic flow, (v(r)). fn the steady 
state, the time-averaged velocity of each individual seg- 
ment of the grafted chains must vanish because they are 
irreversibly grafted to the substrate. This property implies 
J dx (f>h{x)vh.y{x) = 0, which is obeyed in our simulations 
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as shown in Fig. 2a. The drag force, exerted by the x- 
dependent, hydrodynamic flow field, couples the 

motion in the parallel and perpendicular directions and 
gives rise to a correlation between the position, x, of a 
brush segment and its velocity, Vy. A brush segment that 
is farther away from the grafting surface experiences a 
larger drag force and has, on average, a larger parallel 
velocity, Vh,y, than the same segment when it is located 
closer to the grafting surface and exposed to a smaller 
(inverted) drag. This effect is illustrated in Fig. 3. In 
Fig. 2b, one clearly observes that the very few segments 
of the brush that extend far into the melt, x/Rc > 2, 
are dragged along with the flow and have a large posi- 
tive velocity {v\y^y{x) > ^{vy){x)). Since they are grafted, 
they must compensate this by a negative velocity when 
they are closer to the grafting surfaces. This mechanism 
results in the cyclic motion of the individual molecules. 
Our simulations show that this cyclic motion is not nec- 
essarily a consequence of the non-linear elasticity of the 
chain molecules [14] and that it also persists in moder- 
ately dense brushes. Averaging the velocity of all brush 
segments in the vicinity of the surface, we obtain a neg- 
ative velocity and, since the density of the melt chains is 
low at the surface, the negative velocity inside the brush is 
sufficient to invert the total, mass averaged velocity profile 
(v) (cf. Fig. 2a). 

For the choice of parameters, the velocity profile of the 
free chains of the melt does not exhibit an inversion of the 
flow direction at the surface. Since the free chains only 
weakly penetrate the brush, the larger part of the extended 
molecule is exposed to flow in the positive direction at 
the centre even if a single segment is located close to the 
surface. In this case, the drag force on this near-surface 
segment is compensated by the bonding force (as sketched 
in the middle of Fig. 3) and no inversion of the flow is 
observed. To observe the inversion of the flow direction 
of freely moving particles at the surface, we have added 
1% of monomeric tracer particles (i.e., 6554 beads), which 
have identical size than the polymer segments and friction 
coefficient, ( = 1. In Fig. 2a, we demonstrate that the 
velocity profile of these free tracer particles exhibits flow 
inversion at the surface and that such particles potentially 
can be used to experimentally detect the inversion of the 
flow direction at the surface. 

We expect that the inversion of the velocity is charac- 
teristic for intermediate grafting densities. At very small 
grafting densities, the flow at the surface is dominated 
by the velocity of the solvent and the cyclic motion of the 
grafted chains will not be sufficient to invert (v) . Contrary, 
if the brush was strongly stretched, the width of the brush- 
melt interface would be narrow, the height fluctuations of 
brush segments would be small and - in the Lennard- Jones 
model and experiments - topological constraints would be- 
come more important. These effects would shift the in- 
version zone away from the grafting surface and tend to 
reduce the inversion of the flow. 

To verify that the inversion of velocity profile at the 




Figure 4: Mass averaged, total velocity profile, (vy), and ve- 
locity of brush and melt segments for Couette flow. Panel (a) 
highlights the total velocity at the surface and shows details of 
(l>hVh,y for a system that additionally contains 1% tracer parti- 
cles (cf. Fig. 2) . Panel (b) shows the velocity profiles across half 
of the film. The dashed line is a fit to macroscopic predictions 
(cf. Eqs. (3)). The arrow marks the positions, xc, at which the 
extrapolation of the macroscopic hydrodynamic velocity profile 
vanishes. 

brush-coated surface is independent from the way the flow 
is generated we have studied Couette flow, moving the 
grafted chain ends at the top surface with constant velocity 
'jv^Dx- The velocity profile across the pore is shown in 
Fig. 4. ^ The linear behaviour at the centre corresponds 
to macroscopic hydrodynamics 

fhydro,a(a;) = 7 (a; - xq) (3) 

but at the surface the flow direction is inverted, xq de- 
notes the location where this macroscopic linear velocity 
proflle vanishes. If we used this as the definition of an 
effective slip length, the slip length would be negative [9]. 

To validate this observation, we additionally perform 
MD simulations of a standard Lennard- Jones bead-spring 
model. Fig. 5a presents the density profiles of brush and 
melt segments while panel b presents the total velocity 
proflles for Poiseuille and Couette flow. The inversion of 
the flow direction inside the brush is confirmed. Impor- 
tantly, the MD simulations do not invoke any approxima- 
tion: They do not ignore inertia effects and do not re- 
place intermolecular interactions by quasi- instantaneous 
external fields [26]. They capture the local packing of a 
dense polymeric fluid (see density oscillations in Fig. 5a), 
the flnite extensibility of bonds and, importantly, the non- 
crossability of the chain molecules. 

^Only the bottom half of the slit is shown in Fig. 4 Due to Galilean 
invariance, the flow velocity at the top wall is positive and segments 
at the top wall move in the positive direction with a higher velocity 
than the top wall. 
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Figure 5: a) Density profile of the brush-melt interface ob- 
tained from the MD simulations using a bead-spring model 
{N = 32, Re = e.SlcTLj, aRl = 8.9, pocr£j = 0.69, V^iv" = 
5.95), fcsT/eLj = 1.68). b) Mass averaged velocity profile (vy) 
for Poiseuille flow with -FcxCtlj = 0.025eL,7 and Couette flow 
with Wi= 46. The data for Couette flow have been scaled by a 
factor 0.4. The arrows mark xp and xc- The inset highlights 
the velocity at the surface. 

Utilising the simulation results of Poiseuille and Cou- 
ette flow wc try to extract the parameters of the hydro- 
dynamic boundary condition [23] that serves to match 
the detailed microscopic velocity profile in the vicinity 
of the brush-coated surface with macroscopic hydrody- 
namics, Vhydro- The Navier slip condition [21] assumes 
that the viscous stress a^y'^ = Ty '^"''^^""'" calculated from 
the macroscopic hydrodynamic velocity profile, Vhydio, 
equals the friction stress — ^Whydro,^ at a distance, 
Xh- Sh denotes the slip length and xy, the effective po- 
sition of the hydrodynamic boundary. Both parameters 
are assumed to be materials properties that characterise 
the coated surface independently from the strength and 
type of flow. From Eqs. (2) and (3) one obtains 6h — 
■y/ (xp ~ xc){D — xq — Xp) and Xh — xc+Sh- Thus, a nec- 
essary condition for simultaneously describing both types 
of flows by the same parameters, Xb and (5b, is xp > xc 
which is neither fulfilled by the profiles extracted from the 
SCMF simulations with — lORe (see Figs. 2 and 4) nor 
by those obtained from MD simulations of the bead-spring 
model with = G.URe (cf. Fig. 5). 

Finally, in Fig. 6 we analyse the fluctuations of the an- 
gle, (fi, between the projection of the chain's end-to-end 
vector on the grafting surface and the flow direction (see 
inset) for different shear rates. Similar to the behaviour 
of isolated chains [18,20], the wings of the distribution 
for moderate Weissenberg numbers are compatible with 
P{f) ^ l/sin^(/5. This non-Gaussianity of the distri- 
bution is corroborated by the MD results. At compa- 
rable Weissenberg number, the effect is smaller for the 
Lennard- Jones brush because of the higher grafting den- 
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Figure 6: Probability distribution of the angle, ip (see lower 
sketch). For moderate shear, the distribution is sharply peaked 
around </3 = and its wings are compatible with P{}f) ~ 
1 / sin^ (ip) (dashed line) . SCMF simulation data are obtained 
for D:c = 6-Re and Ly = = 2Re. The time step At for 
the highest shear rate, Wi= 34.5, has been decreased by a fac- 
tor 0.1. The MD simulation data of the bead-spring model 
with Wi= 46 are presented by flUed symbols. The upper in- 
set presents the width of the distribution as a function of the 
Weissenberg number. 

sity and non-crossability. The full width at half max- 
imum A(pFWHM of the distribution is compatible with 
ifSFWHM Wi""'^ in agreement the corresponding pre- 
dictions for isolated chains and low shear rates. 

Conclusion. — In summary, dynamic SCMF and MD 
simulations of polymer brushes in an explicit solvent show 
that the motion of chains in a dense brush is similar to the 
cyclic tumbling of isolated grafted chain in shear flow [17]. 
We find that brush chains are not uniformly tilted by the 
flow but the orientation of the molecular axis with the 
flow direction is characterised by a non-Gaussian distribu- 
tion [18,20]. The individual tumbling motion of the brush 
molecules leads to an inversion of the total, collective flow 
velocity at the grafting surface. It would be interesting to 
quantitatively compare simulations with an implicit sol- 
vent using Brinkman's equation and our simulations with 
an explicit solvent to assess to what extent the brush can 
be described as a static porous medium. The simulation 
also provide microscopic insights into the mechanisms that 
dictate the hydrodynamic boundary condition. The lat- 
ter is important for controlling flow in microfluidic de- 
vices. The cyclic tumbling motion of the individual brush 
molecules could also potentially serve to separate parti- 
cles with different surface affinities by flow or to promote 
mixing close to the grafting surface. 
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